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ABSTRACT 

We report the first detection of the Zeeman effect in the 36 GHz Class I 
CH3OH maser line. The observations were carried out with 13 antennas of the 
EVLA toward the high mass star forming region M8E. Based on our adopted 
Zeeman splitting factor of z = 1.7 Hz mG -1 , we detect a line of sight magnetic 
field of —31.3 ± 3.5 mG and 20.2 ± 3.5 mG to the northwest and southeast of 
the maser line peak respectively. This change in sign over a 1300 AU size scale 
may indicate that the masers are tracing two regions with different fields, or that 
the same field curves across the regions where the masers are being excited. The 
detected fields are not significantly different from the magnetic fields detected 
in the 6.7 GHz Class II CH 3 OH maser line, indicating that CH 3 OH masers may 
trace the large scale magnetic field, or that the magnetic field remains unchanged 
during the early evolution of star forming regions. Given what is known about the 
densities at which 36 GHz CH 3 OH masers are excited, we find that the magnetic 
field is dynamically significant in the star forming region. 

Subject headings: ISM: magnetic fields — masers — polarization — radio lines: 
ISM — stars: formation 



1. INTRODUCTION 



Magnetic fields likely play an important role in the star formation process (e.g.. lTroland et al 



20081 . and references therein). The nature of this role is not yet clear, however, primarily 
due to the scarcity of observational data on magnetic fields in star forming regions. The 
Zeeman effect remains the most direct method for measuring magnetic field strengths (e.g., 



Crutcherlll999l ). Observations of the Zeeman effect in H I and OH thermal lines have revealed 
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the strength of the mag netic field in the lower density envelopes of molecular clouds (e.g., 
Brogan fc Trolandll200ll ). Observations of the Zeeman effect in H 2 masers, on the other 
hand, offer a window into magnetic fields in the highest densit y regions (n ~ 10 9 cm -3 ) of 
star forming regions (ISarma et al.ll2008l ; IVlemmings et al.ll2006l ). 



Interstellar masers, being c ompact and intens e, are extremely effective probes of young 
star forming regions. Recently, IVlemmingsl (120081 ) carried out the first systematic study of 
the Zeeman effect in the 6.7 GHz methanol maser line. Such (6.7 GHz) masers are exam- 
ples of Class II methanol masers, which are known to be probes of the early phases of high 
mass star forming regions. Yet another class of masers , namely Class I m ethanol masers, 
may probe even earlier phases of star forming regions (IPratap et al.l 120081 ). This appears 
to be the case for M8E, which is a high mass star forming region located at I = 6°. 05, 
b = -T.45, at a distance of 1.5 kpc Jsimon et al.lll984l ; iGeoreelin fc Georgelinlll976h . M8E 
is comprise d of a bright infrared source (M8E-IR ), and a radio H II region about 8" to its 
no rthwest (ISimon et al.lll984l ; Ide Wit et al.l 120091 ) . Low spatial resolution CO observations 
bv lWright et al.l (119770 showed that M8E is embedded in a dense molecular core. At higher 
resolution, Mitchell et al.l (119921 ) mapped a bipolar outflow in the 12 CO emission line (i'lsr= 



11 k m s 1 ), with blue a nd red wings spread over 50 km s 1 . Based on their VLTI observa- 



Linz et al.l (120091 ) concluded that M8E-IR contains a central protostar of mass 10-15 



tions, 
M . 

In this paper, we report the first detection of the Zeeman effect in the 36 GHz Class I 
CH 3 OH maser line. The observations were carried out with 13 antennas of the Expanded 
Very Large Array (EVLA) toward the high mass star-forming region M8E, and represent an 
early success of the EVLA. In § [21 we present details of the observations and reduction of 
the data. The analysis involved in extracting magnetic field information from the Zeeman 
effect is given in § [31 The results are presented and discussed in § HI 



2. OBSERVATIONS & DATA REDUCTION 

Observations of the 4_i — 3o methanol maser emission line at 36 GHz were carried out 
using the Expanded Very Large Array (EVLA) of the NRAC0 in two 2 hr sessions on 2009 
July 9 and 25. Thirteen EVLA antennas equipped with the 27—40 GHz (Ka-band) receivers 
were used in these observations. To avoid the aliasing that affects the lower 0.5 MHz of the 
bandwidth for EVLA data correlated with the VLA correlator, the spectral line was centered 



1 The National Radio Astronomy Observatory (NRAO) is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. 



- 3- 



in the second half of the 1.56 MHz wide band. The source 3C286 (J1331+3030) was used 
to set the absolute flux density scale, while the compact source J1733— 1304 was used as an 
amplitude calibrator. 

The editing, calibration, Fourier transformation, deconvolution, and processing of the 
data were carried out using the Astronomical Image Processing System (AIPS) of the NRAO. 
After applying the amplitude gain corrections of J 1733— 1304 on the target source M8E, the 
spectral channel with the strongest maser emission signal was split, then self-calibrated in 
both phase and amplitude in a succession of iterative cycles. The final phase and amplitude 
solutions were then applied on the full spectral-line uv data set, and Stokes / and V image 
cubes were made. In the channel with the maser li ne peak, we recovere d about 70% of the flux 



density reported in the single dish observations of lPratap et al.l (120081 ). Further processing of 



the data, including magnetic field estimates, was done using the MIRIAD software package. 



ANALYSIS 



For cases in which the Zeeman splitting Ai/ Z is much less than the line width Av, the 
magnetic field can be obtai ned from the Stoke s V spectrum, which exhibits a scaled derivative 
of the Stokes I spectrum (IHeiles et al.lll993l ). Here, consistent with AIPS conventions, / = 
(RCP+LCP)/2, and V = (RCP-LCP)/2; RCP is right- and LCP is left-circular polarization 
incident on the antennas, where RCP has the standard radio definition of clockwise rotation 
of the electric vector when viewed along the direction of wave propagation. Since the observed 
V spectrum may also contain a scaled replica of the I spectrum itself, the Zeeman effect can 
be measured by fitting the Stokes V spectra in the least-squares sense to the equation 



V 



T b dl 

al H — 

2 dv 



(ITroland &: Heileslll982l ; ISault et al.lll990l ). The fit parameter a is usually the result of small 
calibration errors in RCP versus LCP, and is expected to be small. In these observations, 
a was of the order of 10~ 4 or less. While eq. (CQ) is strictly t rue only for thermal lines , 
numerical solutions of the equations of radiative transfer (e.g., iNedoluha fc Watsonlll992l ) 
have shown that it gives reasonable values for the magnetic fields in masers also. In eq. (JT]), 
the fit parameter b = zBcosO, where z is the Zeeman splitting factor (Hz mG" 1 ), B is the 



magn etic field, and 9 is the angle of the magnetic field to the line of sight fjCrutcher et al. 



19931 ). The value of the Zeeman splitting factor z for CH 3 OH masers is very sma l l, bec ause 
CH3OH is a non-paramagnetic molecule. Following the treatment of I Vlemmingsl (120081 ) for 
the Zeeman splitting of 6.7 GHz methanol masers, we derive the Zeeman splitting factor 
using the Lande g-factor based on laboratory measurements of 25 GHz methanol masers 
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( jJenlll95ll ). and find z = 1.7 Hz mG -1 . While this will likely introduce a systematic bias 
into the magnetic field results, it remains the best possible estimate pending laboratory 
measurements of the 36 GHz maser line. 



4. RESULTS & DISCUSSION 

In order to demonstrate our detection of the Zeeman effect in the 36 GHz Class I CH 3 OH 
maser line toward M8E, we display the Stokes I and V profiles toward two positions in this 
source in Fig. [1} the two positions are to the northwest and southeast of the maser line peak. 
As described in § [31 we determined magnetic fields by fitting the Stokes V spectra in the 
least-squares sense using equation (JTJ). The values of the fit parameter b (see eq. [1]) obtained 
from this fit are b = —53.2 ± 6.0 Hz for the northwest position, and b = +34.4 ± 5.9 Hz for 
the southeast position. Using the Zeeman splitting factor 2 = 1.7 Hz mG^ 1 discussed in § 
[3] above, this gives B\ os = —31.3 ± 3.5 mG for the northwest position, and 20.2 ± 3.5 mG 
for the southeast position. By convention, a negative value for £? los indicates a field pointing 
toward the observer. 

The line-of-sight magnetic field has opposite signs at the two positions for which Stokes 
I and V profiles are shown in Fig.[U with values of —31.3 mG and 20.2 mG, respectively. The 
observed change in the sign of B\ os at these two positions, together with a slight asymmetry 
in the maser line profiles at each position, indicates that we are observing at least two masers 
that are very close in position and velocity. The masers are marginally resolved in our C- 
configuration observations, otherwise the opposite magnetic fields would sum to zero. EVLA 
B- or A-configuration observations will be necessary to fully resolve the maser components. 
Such observations will not be possible until late 2010 or early 2011. The observed change 
in sign of B\ os occurs over a size scale of 0".9, equal to 1300 AU (assuming the distance to 
M8E is 1.5 kpc). This may mean that the clumps where the 36 GHz maser is being excited 
come from two different regions where the field is truly different. Alternatively, it may mean 
that the field lines curve across the region in which the masers are being excited, so that 
the line-of-sight field traced by one maser is pointed toward us, whereas that traced by the 
other maser is pointed away from us. 

Since this is the first detection of the Zeeman effect in the 36 GHz CH 3 OH maser line, 
we will address the issue of beam squint, a topic which seems to arise every time a Zeeman 
detection is mentioned. Beam squint arises because the RCP and LCP beams have slightly 
different pointing centers. Therefore, if a velocity gradient is present across an extended 
source, beam squint can give rise to a false Zeeman pattern. Since masers are compact 
sources, beam squint should not be an issue here. Moreover, since the signs of the detected 
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field are reversed at the two positions, it would be difficult to manufacture a velocity gradient 
that would reverse the signs of the fields as observed. 

Our detected B\ os values in the 36 GHz Class I CH 3 OH maser line toward M8E are 
similar to the val ues obtained in the 6. 7 GHz Class HCH 3 OH maser line observed by 
Vlemmingsl (120081 ) toward other sources. IVlemmingsl (120081 ) detected significant magnetic 
fields in 6.7 GHz methanol masers with the 100 m Effelsberg telescope in 17 sources, with 
an average value of 23 mG. Since Class I and Class II CH 3 OH masers trace different spatial 
regions, this may indicate that the magnetic field is the same over the extent of these two 
regions, giving rise to the possibility that CH 3 OH masers trace the large-scale magnetic 
field in the star forming region. Alternatively, if Class I masers occur in the very early 
stages of star formation (before the formation of an ultracompact H II region), and Class II 
masers occur later on, the similarity in B\ os for these two classes, if true, may indicate that 
the magnetic field strength remains the same during the early stages of the star formation 
process. Certainly, more observations are necessary, especially in regions where both Class I 
and Class II masers are known to occur, in order to address the possibilities noted above. 

Finally, we use our observed value for the magnetic field to compare the magnetic 
and dynamical energies in these masing regions. The magnetic energy de nsity is given by 
3Bf os /8n. Here, we have used the relation B 2 = 3Bf os that was obtained by ICrutcherl (119991 ) 
on statistical grounds. For B\ os 25 mG from these observations, the magnetic energy 
density will be 8 x 10 -5 ergs cm -3 . The kinetic energy density (thermal and turbulent) is 
given by (3/2)mna 2 , where a = Au/(81n2) 1 / 2 is the velocity dispersion, the mass m = 2.8m p 
(assuming 10% He; m p is the proton mass). Since masers occur only in special directions 
along which they have developed the required velocity coherence, the velocity dispersions in 
the masing region may be g reater than th at traced by masers. If we use Av = 5 km s _1 - 



McKedll999l ). we get a kinetic ene r gy de nsity equal to 3 x 10 



a reasonable value (see, e.g., 
ergs cm -3 for n = 10 5 cm -3 , the density quoted in lPratap et al.l (120081 ) at which maser action 
in the 36 GHz line is maximized. Moreover, suppose we find fields of similar strength in the 
44 GHz CH3OH maser line in a future observation, 44 GHz being the other prominent Class 
I CH3OH maser line. Then, c onsidering the dens ity n = 10 6 cm -3 at which maser action in 



the 44 GHz line is maximized (IPratap et al.ll2008l ). we will get a kinetic energy density equal 
to 3 x 10~ 7 ergs cm -3 . Either way, the magnetic energy is dominant by at least two orders 
of magnitude, indicating that the magnetic fields must be dynamically significant in these 
regions. 
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5. CONCLUSIONS 

We have detected for the first time the Zeeman effect in the 36 GHz Class I CH3OH 
maser line toward the high mass star forming region M8E. The detected values for the line-of- 
sight magnetic field are —31.3 ±3.5 mG and 20.2 ±3.5 mG. There may be systematic bias in 
these values due to the assumed Zeeman splitting factor of z — 1.7 Hz mG" 1 , which is based 
on laboratory experiments on the 25 GHz CH3OH maser line. Our detect ed values are simila r 



to the magnetic fields detected in 6.7 GHz Class II CH3OH masers by IVlemmingsl (120081 ). 
This suggests that CH3OH masers may trace the large-scale magnetic field; alternatively, the 
magnetic field may remain the same during the early evolution of the star forming region 
between the time when Class I and Class II masers are excited. Our detected values for Bi os 
imply that the magnetic field is dynamically significant to these regions. 

Our detection gives rise to a host of questions that need further observational and 
theoretical work. We are planning future observations of the Zeeman effect in more 36 
GHz masers and at higher resolution. We are also planning for observations of 44 GHz 
masers with the aim of detecting the Zeeman effect. As mentioned above, the 44 GHz line 
constitutes the other prominent Class I CH3OH maser line. We are also planning to carry 
out high angular resolution observations of the Zeeman effect in regions containing both 
Class I and Class II masers, in order to compare the strength of magnetic fields between the 
regions traced by these classes. From theorists, we seek detailed models for Class I masers in 
order to understand their properties better, especially the range of densities at which they 
are excited. Also, there is a pressing need for laboratory studies of the Zeeman splitting 
coefficient for CH 3 OH masers. The resolution of some or all of these issues will open up yet 
another important window into the early stages of high mass star forming regions. 
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Table 1. Parameters for EVLA Observations 



Parameter 


Value 


Observation Dates 


2009 July 9 & 25 


Configuration 


C 


R.A. of field center (J2000) 


18 h 04 m 53 s .3 


Decl. of field center (J2000) 


-24°26'42".0 


Total Bandwidth 


1.56 MHz 


No. of channels 


256 


Channel Spacing 


0.051 km s" 1 


Time on Source 


2 hr 


Rest Frequency 


36.16929 GHz 


Velocity at band center 3. 


13.7 km s _1 


Target source velocity 


11.2 km s" 1 


Hanning Smoothing 


No 


FWHM of synthesized beam 


1.76" x 0.58" 




PA. = -7.80° 


Line rms noise b 


18 mJy beam -1 



a The line was centered in the second half of the 
1.56 MHz band in order to avoid aliasing (see §|2j). 

b The line rms noise was measured from the 
stokes I image cube using maser line free channels. 
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V LSR ( km S ^ V LSR ( km S l ) 

Fig. 1. — Stokes / (top— histogram) and V (bottom— histogram) profiles of the maser toward 
the (a) northwest and (b) southeast of the maser line peak. The curve superposed on V in 
each of the lower frames is the derivative of I scaled by a value of B\ os = —31.3 ± 3.5 mG in 
(a), and B los = 20.2 ± 3.5 mG in (b). 



